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ABSTRACT 


One of the key questions in solar physics that remains to 
be answered concerns the strength and the distribution of 
the magnetic fields at the base of the convection zone. 
The flux tube dynamics requires that toroidal fields of 
strength as large as 100 kilogauss be present at the base of 
the convection zone. The kinetic-magnetic equipartition 
argument leads to smaller field strengths. For possible 
detection of these relatively small (compared to pressure 
effects) fields by helioseismic methods it is important to 
know the range of the field strengths and their distribu¬ 
tion. 

We estimate a range for the toroidal magnetic field 
strengths at the base of the convection zone using dy¬ 
namo simulations in a spherical shell. These simulations 
involve the distribution of rotation provided by helioseis¬ 
mic inversions of the GONG and MDI data. Combining 
the simulations with the observed line-of-sight surface 
poloidal field we extract the spatial pattern and magni¬ 
tude of the mean toroidal magnetic field at the base of the 
convection zone. 


1. INTRODUCTION 


The inner distribution of the solar magnetic field is poorly 
known. Despite the striking regularity of the solar activ¬ 
ity, we still lack an understanding of its causes and the 
mechanisms by which it operates. 

The layer that we are able to directly observe the mag¬ 
netic fields on is the photosphere. From observations of 
sunspots we know that magnetic fields with intensities 
of the order of thousands of Gauss are produced some¬ 
where below the photosphere and can live for months be¬ 
fore decaying away. The topological and magnetic struc¬ 
ture of sunspots (as expressed by Hale’s and Joy’s laws) 
suggest that they harbor intense, toroidal magnetic flux 
tubes. The poloidal component of the solar magnetic field 
is much weaker compared with the toroidal field. 

The intensity and spatio-temporal distribution of mag¬ 
netic fields inside the Sun represent a challenge for any 
theory explaining the origin of solar activity. On one 
hand, the highly localised sunspot field points toward 


local dynamics. On the other, the much more diffuse 
poloidal field more naturally fits into a mean-field model. 
One needs to combine two approaches to explain the so¬ 
lar cycle: the flux tube dynamics and the mean field 
dynamo. Thin flux tube dynamics has been very suc¬ 
cessful in explaining features connected to sunspots, see 
e.g Choudhuri & Gilman 1987, Fan Fisher & De Luca 
1993, Caligari, Moreno-insertis, & Schiissler 1995. The 
periodic migration of sunspots toward the equator and 
the field reversal cannot be explained by this dynam¬ 
ics alone. The explanation of the reversals of the field 
and of the equator-ward migration of the activity were 
in fact one of the great achievements of the mean field 
dynamo theory originated by Parker, Steenbeck, Krause 
and Radler (Parker 1955, Steenbeck, Krause and Radler 
1969). Yoshimura (1975a) proved that waves mainly 
propagate along the isosurfaces of angular velocity. With 
the present knowledge of solar differential rotation, this 
would imply that inside the bulk of the convection zone 
waves travel radially outwards, not giving rise to equator- 
ward migration. Migration may be restored in the shear 
layer at the base of the convection zone. 

While it is possible to generate high magnetic fluxes by 
means of differential rotation, to store them for enough 
time to allow for their intensity to build up would not be 
feasible in the convection zone because such magnetic 
flux tubes would erupt in a timescale of months (Parker 
1975, Moreno-insertis 1986). Helioseismology tells us 
that around 0.7 R 0 , a sharp radial change in the solar 
rotation curve happens in a layer whose thickness can 
be as small as 0.02 R 0 (see Christensen-Dalsgaard et al. 

1991, Basu & Antia 1997, Kosovichev 1996, Corbard et 
al. 1998, Charbonneau et al. 1999.). 

The same layer may allow for fields up to 10 5 Gauss to be 
stored, what would be needed for thin flux tube dynam¬ 
ics to work in the case of the Sun (Moreno-insertis et al. 

1992, Ferriz-Mas & Schiissler 1994) 

Following Ivanova & Ruzmaikin (1977), Parker (1993) 
discussed a model of mean field dynamo that includes 
a sharp gradient of turbulent diffusivity and two distinct 
location for the sources of the magnetic field, that is dif¬ 
ferential rotation and helical turbulent motions. Stronger 
toroidal fields can in fact be produced in the region where 
diffusivity is smaller, just below the convection zone. 
Separating the shear layer and the source of the alpha- 
effect, moreover, would allow for alpha not be quenched 
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by the strong underlying magnetic field. A very thorough 
investigation of this kind of models has been carried out 
by Charbonneau and MacGregor, 1997. We will use in 
our Profile II a similar setup to one of those discussed in 
the aforementioned paper. 

In the following, we are going to study how different pro¬ 
files of the a-effect may influence the spatial distribution 
of fields. We shall consider the profile of diffusivity and 
the rotation curve as given. 

Mean field, kinematic dynamo cannot predict the abso¬ 
lute values of the generated magnetic fields, only their 
ratio. The measured mean radial field at the surface can 
then be used to infer the value of the toroidal magnetic 
field in the interior, once the ratio is known. 


2. DYNAMO MODEL 


We consider axisymmetric solutions of the mean field dy¬ 
namo. We assume for the solar rotation a simple analyt¬ 
ical fit to the profile reconstructed by helioseismic where 
the surface rotation curve 

O s = 0£q (l + ai cos 2 0 + 02 cos 4 6 ) (1) 

is made smoothly match the core rotation in a layer 
of thickness 0.2Rq at O.692R 0 . Equatorial rotation is 
Ofiq = 2.865 x 10~6s _1 and the core rotation taken 
as the value of O s at 30° latitude. 9 here is colatitude. 
a\ = —.126 and 02 = —0.159. The radial profile of 
dLl/dr at the equator is shown in solid line in Figure [lj. 
The sign of this gradient is opposite at higher and lower 
latitudes. Turbulent diffusivity p is constant throughout 
the Convection Zone and we have it drop a factor 10 
to 50, in a layer of thickness 0.2R o at a location ei¬ 
ther coincident with that of the rotational shear layer or 
slightly above it, at O.713R 0 . Both these values have 
been worked out in the context of Helioseismology, see 
Christensen-Dalsgaard et al. 1991, Basu & Antia 1997, 
Kosovichev 1996, Corbard et al. 1998, Charbonneau et 
al. 1999. In reality the drop in turbulent diffusivity is 
estimated to be of the order of 10 6 . In Figure [I] the diffu¬ 
sivity profile along the radial direction is represented by 
the solid dashed line, its gradient marking the bottom of 
the Convection zone. In solid line, the radial derivative of 
the rotation which defines the shear layer is plotted. 

The mean field dynamo equation, written in terms of 
scalar potentials (Krause & Radler, 1980) is solved nu¬ 
merically in a r — 6 meridional semi-disk. Second order 
finite differences in space and a third order Runge-Kutta 
scheme for time advance are used. In most of the runs a 
grid of 60 x 40 is used. Ideal conductor boundary con¬ 
ditions are used at the interface with the core and radial 
field condition is implemented at the surface. Regualrity 
condition are imposed on the axis. 

No model for non-linear a-quenching has been used. 
This is consistent with our assumption that the flow field 
is not influenced by the magnetic field. 


2.1. The a-effect 


Given the lack of constraints on the form of the a-effect, 
we shall discuss how different alphas influence the spatial 
distribution of the the toroidal and poloidal fields. Fig¬ 
ure [I] shows the distribution in radius of a. In all models, 
a latitudinal dependence of cos (6) is included. This is 
standard in dynamo theory and reflects the property of 
the alpha-effect of being antisymmetric with respect to 
the Equator. In case IV, an additional factor of sin 2 (9) 
is present. This has been discussed by Rudiger & Bran¬ 
denburg, 1995, and has the effect of shifting the magnetic 
field patterns closer to the equator, which better repro¬ 
duces the observed patterns of sunspot migration. 


Profile I, the diamonds in Figure [l], is maximum in the 
bulk of the Convection Zone, at 0.82 Rq in our case. It 
is null at the surface. This form of the a-effect takes into 
account the influence o the rotation on helical turbulence, 
see Zeldovic, V. Ruzmaikin, Sokoloff 1983. Profile II, 
squares, is sharply peaked at 0.7 f? 0 , just above the shear 
layer, and represents an Interface Dynamo model. As al¬ 
ready mentioned, this profile was used by Charbonneau 
and MacGregor (1997). Profile III is, instead, non van¬ 
ishing only in the outer shell of the Convection Zone. 
This kind of “surface dynamo” has been discussed in con¬ 
junction with meridional circulation, see e.g. Choudhouri 
et al. 1995. Profile IV, triangles, has a constant value 
throughout the whole Convection Zone. It drops to zero 
close to the surface and at the base of Convection Zone. 
This profile has often be adopted in the literature. The last 
case considered, V, a has two contributions: one com¬ 
ing from the bottom of the convection zone where it is 
negative and proportional to the gradient of p, combin¬ 
ing the effects due to the the decrease in the intensity of 
the turbulence and the stratification. The other, coming 
from the convection zone, is taken to be the same as in I. 
Yoshimura (1975b) has used this form for the a-effect. 


a(r) 

r a 

ao 

7 

T 

rmax 

Bt/Br 

i 

.75 

4.1 

8.2 

.026 

.71, .84 

310 

ii 

.7 

27.5 

19 

.047 

.67, .80 

950 

hi 

> .8 

4.1 

18 

.008 

.78, .92 

36 

IV 

.7 -7- .98 

4.23 

5 

.015 

.9 

73 

V 

.75 

5 

180 

.020 

.71 

360 


Table 1. Location and intensity of the maximum toroidal 
magnetic field for different choices of the profile of the 
a-effect. The location where the profile is peaked r a 
is given, along with the intensity ao- 7 and T are the 
growth rate and the period of the solutions. B r is the 
maximum value of the radial photospheric field averaged 
over time. Bf is the maximum value of the time average 
of the toroidal magnetic field intensity. When two values 
are given, the first one represents a secondary maximum 
close to the shear layer, see Figure rmax A the location 
where the maximum is located, in units of Rq. Numbers 
are adimensional. Lengths have been scaled to Rq and 
time to Rq/t], where p is the value of turbulent diffusivity 
in the convection zone. 
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The a —effect, diffusivity and rotation profiles. 



Figure 1. Radial distribution of the a-effect for the cases I-V. Values are not to scale. See Table 0/or the actual values 
used in the present simulations. The solid lines, dashed and continuous, represent the radial profile of the turbulent 
magnetic diffusivity p and the radicd gradient of the rotation rate dfl/dr at the equator, respectively. They define for our 
model the Convection Zone and the shear layer ( the Tachocline). 




Figure 2. Surface plot of toroidal and radial field intensity integrated over time for Profile I. Values are scaled to the 
maximum of the radial surface field. 


3. DISCUSSION OF THE RESULTS 


Except for the case of Profile III, where a is concen¬ 
trated in the surface layer, all the case studied display a 
local maximum of the toroidal field in the shear layer. 
Except for case V, this is not an absolute maximum 
which is instead achieved in the convection zone, within 
0.8 -r 0.9 Rq, see Figure 0 and Figure |[ 

The ratio of the maximum toroidal field to the surface ra¬ 
dial field can vary from a few tens as in the case of the 
surface a-effect, case III, to a thousand, as in the case of 


the interface profile II. Both the profile I and V have ra¬ 
tios of the order of a few hundreds. No symmetry across 
the equator has been imposed on the solution which dis¬ 
play a North-South asymmetry. Except for case of the 
surface a, the symmetry of the solution is mainly dipolar. 
Considering a mean radial surface field of the order of a 
Gauss (Schlichenmaier and Stix 1995) the range of ratios 
that we found would lead to an estimate of the large-scale 
mean toroidal field of the order of 10 3 Gauss in the most 
favorable cases I,II,V. Both the surface (III) and top-hat 
(IV) profiles have smaller ratios of the toroidal to the sur¬ 
face fields. 
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Figure 3. Toroidal fields for Profile II-V. Values are scaled to the maximum of the radial surface field. 


If one assumes that solar activity is originated in the shear 
layer, then it is possible to draw a butterfly diagram with 
the temporal evolution of the toroidal field near the shear 
layer. All those models show butterfly diagrams with ac¬ 
tivity at higher latitude then observed. This is a known 
feature of many dynamo models. 

We have shown that different assumptions about the a- 
effect give rise to very different spatial distributions of the 
magnetic fields inside the Sun. This messsage could also 
be read in reverse: should we be able to probe the field 
deeply inside the Sun, we could have information about 
the nature of the regeneration mechanism expressed by 
the a-effect. 
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